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ABSTRACT 


1.0. Introduction 


A time-lapse seismic survey compares two or more seismic surveys acquired at different stages of hydrocarbon 
production (Stammeeijer and Hatchell, 2014). Time-lapse seismic monitoring of production-induced changes in a 
reservoir and the surrounding rocks over time has the basic aim of mapping reservoir compartments and 
subsurface rock deformation, monitoring fluid movement and pore pressure changes, identifying by-passed oil and 
planning for future production performance, and this has been challenging (Landrg, 2001; Trani et al., 2011; 
MacBeth et al., 2018). Such integrated geological, geophysical and petroleum engineering methods are needed to 
create a well-resolved picture of a producing reservoir. Pore pressure depletion creates changes in the stress and 
strain fields of the rock material both inside and outside the reservoir (Prioul et al., 2004; Fuck et al. (2011)). In 
addition to leakage of hydrocarbons, hazards are associated with wells crossing reactivated faults (Alsos eg al. 


2002). Reactivated faults can also have close relationships with earthquakes (Zoback and Zincke, 2002). 


Many literatures have reported interesting seismic timelapse studies in the oil prolific area of the Niger Delta 
Basin, Nigeria. Furthermost, these studies majored on timelapse feasibility studies and reservoir monitoring and 


management (Aniwetalu et al., 2017; Uko and Otugo, 2016). 


Moreover, literatures in public domain are very scarce on seismic timelapse, in Nigeria, with the main objective of 
determining effect of hydrocarbon production on subsurface structures and faults reactivation (Igwenagu et al, 
2021). Fuck et al, (2011) introduced an analytic geomechanical model with third-order elasticity to predict 


travel-time shifts due to stress-induced seismic velocity heterogeneity and anisotropy within and surrounding a 


ISSN: 2581-5059 OPEN@ACCESS 116 


S Mediterranean Journal of Basic and Applied Sciences (MJBAS) 
Volume 7, Issue 2, Pages 116-125, April-June 2023 


PROMOTING RESEARCH & KNOWLEDGE 


compacting reservoir. Angus et al, (2015) extended and calibrated a micro-structural rock-physics approach to 
model the nonlinear elastic response due to varying effective stress and hence predicted induced changes in 


seismic velocity and anisotropy. 


This study is taken from Fuba Field, Depobelt, Niger Delta, Nigeria. The ultimate deliverable of this study was the 
subsurface structures in the area, using 3-D seismic time-lapse data. The major components of this study are: (a) 
Well Correlation performed in order to determine the continuity of the reservoir sand across the field. (b) Seismic 
Interpretation which involves well-to-seismic tie, fault mapping, horizon mapping, time surface generation and 
depth conversion. This aids in giving more insight into the subsurface structures in onshore Fuba field Niger-Delta, 


Nigeria, using 3-D seismic time-lapse data. 
== 2.0. Location and Geology of the Study Area 


The proposed study area Fuba Field is located in the onshore Niger Delta region. Figure 1 shows the map of the 
Niger Delta region showing the study area. The Niger Delta lies between latitudes 4° N and 6° N and longitudes 3° 
E and 9° E (Whiteman, 1982). The Delta ranks as one of the major oil and gas provinces globally, with an 
estimated ultimate recovery of 40 billion barrels of oil and 40 trillion cubic feet of gas (Adegoke et al., 2017). The 


coastal sedimentary basin of Nigeria has been the scene of three depositional cycles (Short and Stauble, 1967). 


Figure 1. Map of Niger Delta Oilfields showing the location of Fuba Field 


The first began with a marine incursion in the middle Cretaceous and was terminated by a mild folding phase in 
Santonian time. The second included the growth of a proto-Niger delta during the Late Cretaceous and ended in a 
major Paleocene marine transgression. The third cycle, from Eocene to Recent, marked the continuous growth of 
the main Niger delta. A new threefold lithostratigraphic subdivision is introduced for the Niger delta subsurface, 
comprising an upper sandy Benin Formation, an intervening unit of alternating sandstone and shale named the 
Agbada Formation, and a lower shaly Akata Formation. These three units extend across the whole delta and each 
ranges in age from early Tertiary to Recent. They are related to the present outcrops and environments of 
deposition. A separate member of the Benin Formation is recognized in the Port Harcourt area. It is Miocene- 
Recent in age with a minimum thickness of more than 6,000ft (1829 m) and made up of continental sands and 


sandstones (>90%) with few shale intercalations (Horsfall et al., 2017). Subsurface structures are described as 
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resulting from movement under the influence of gravity and their distribution is related to growth stages of the 
delta (Ochoma et al., 2020). Rollover anticlines in front of growth faults form the main objectives of oil 


exploration, the hydrocarbons being found in sandstone reservoirs of the Agbada Formation. 
“= 3.0. Materials and Methods 
3.1. Normalized Root Mean Square (NRMS) Repeatability (RRR) 


Pre-stack time migrated full-offset 3D and 4D stacks were available. The success of time-lapse reservoir 
monitoring depends on removing the non-repeatable effects such as configurations, seasonal changes, atmospheric 
temperatures, tides, elastic properties of the overburden, compaction, multiples, and rock heterogeneities (Vedanti 
et al., 2009; Varela et al., 2006). Obstructions, weather patterns, cost constraints, and maritime traffic can also 
influence the survey orientation. Having considered the above sources of error in repeatability, the acquisition 
system itself, positioning accuracy, receiver sensitivity/calibration, and source calibration must also be looked 
into. All these sources of error were handled through normalized root mean square (NRMS) analysis. The NRMS 
value is simply the RMS amplitude of the difference, normalized by the average of the RMS amplitudes of the 
baseline data Baseline and monitor data Monitor (Kragh and Christie, 2002): 


2rms (Monitor —baseline ) 


NRMS = (1) 
rms(Monitor )+ (Baseline ) 
2 
>: a 
rms = 
- (2) 
The summation is over N number of all samples x;(i = 1, 2, - - -, N) in the time window. 


The NRMS value is a measure of non-repeatability. If NRMS = 0, the data are perfectly repeatable. Typical “good” 
values of NRMS quoted in the literature range from 0.1 to 0.3 [10% to 30% non-repeatability] (Johnston, 2013). 


Monitor 4D seismic volume was subtracted from the Base 3D volume, and the differences was interpreted to 
determine the areas of the field that have been changed during production. Areas of the areas in the field where 


there have been changes were analyzed and compared to production activity in those areas. 
3.2. Well-to-Seismic Ties 


Well correlation is the first stage of the pre-interpretation process. The process of well correlation involves 
lithologic description, picking top and base of sand-bodies, fluid discrimination and then linking these properties 
from one well to another based on similarity in trends. In between these two lithologies in the subsurface, the 
gamma ray log is often used. Correlation of reservoir sands was achieved using the top and base of reservoir sands 
picked. The correlation process was possible based on similarity in the behavior of the gamma ray log the Niger 
Delta; the predominant lithologies are sands and shales. In order to discriminate shapes. Also, the thickness of the 
shale bodies overlying and underlying the sand body is considered during Correlation. After defining the 


lithologies, the resistivity log was used for discriminating the type of fluid occurring within the pores in the rocks. 


There are six basic steps involved in seismic interpretation relevant to this study and they include; Well-to-seismic 


ties, Fault Mapping, Horizon mapping, Time surface generation, Velocity Modelling and Depth Conversion. 
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Well-to-seismic tie is a process that enables the visualization of well information on seismic data. For this process 
to be achieved, the following are basic requirements; checkshot, sonic log, density log and a wavelet. The sonic 
log, which is the reciprocal of velocity, was calibrated using the checkshot data. The calibration process is 
necessary in order to improve the quality of the sonic log because the sonic log is prone to washouts and other 
wellbore related issues. The results of calibrating the sonic log with the checkshot gives a new log called the 


calibrated sonic log. 


The calibrated sonic log is used along with the density log to generate an acoustic impedance (AI) log. The 
acoustic impedance log is calculated for each layer of rock. The next step involves generating the reflectivity 
coefficient (RC) log. The RC is calculated and generated using the AI log. The RC log generated is then convolved 
with a wavelet to generate a synthetic seismogram which is comparable with the seismic data. The statistical 
wavelet utilized for convolution is extracted from the seismic data. The synthetic seismogram was generated for 
every well that had checkshot, density and sonic log. The reflections on the synthetic seismogram were matched 
with the reflections on seismic data. The mathematical expressions that governs the entire well-to-Seismic tie 


workflow are presented below; 


Al = pv (3) 
RC = P2V2—PiVy4 (4) 
P2V2t p14 
: : P2V2—Pily 
Synthetic Seismogram = ————— * wavelet 5 
8 P2V2tP1Vy (6) 


Where p= Density, v =Velocity, AI= Acoustic impedance and RC = Reflection coefficient 


Faults were identified as discontinuities or breaks in the seismic reflections. Faults were mapped on both inline and 
cross-line directions. Horizons are continuous lateral reflection events that are truncated by fault lines. The horizon 
interpretation process was conducted along both inline and crossline direction. At the end of the horizon mapping, 
a seed grid is generated which serves as an input for time surface generation. Time surfaces were generated using 


the seed grids gotten from the horizon mapping process. 
3.3. Velocity Modelling 


Three velocity models were generated in this study and utilized for depth conversion. The velocity models 
generated includes Linear velocity function (average velocity), second order polynomial and third order 
polynomial velocity models. The linear velocity function, second order polynomial and third order polynomial 
velocity models generated were used separately to depth convert the time surfaces of the reservoirs of interest. 
Uncertainty known as residuals were estimated at well points for the various velocity models. The model with the 


least residual was preferred as most suitable for converting surfaces from time to depth. 

“= 4.0. Results and Discussion 

4.1. Production Data 

The Production and Reservoir Pressure Reports are presented in Figures 2 and 3. Production decline of 1192.21 


bbl/day resulted from pressure decline of 95.50bar/year or 0.062psi/year. 
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4.2. Normalized Root Mean Square (NRMS) Repeatability (RRR) and 4D Response 


The NRMS of 0.38 has been achieved, in this study, implying very good repeatability when considering the quality 
of data and acquisition difference. The seismic time-lapse difference between the base and monitor surveys was 
successfully extracted. The fact of the difference from the monitor implies the existence of production induced 


effect and acquisition, environmental and processing noises, hence the 4D or time-lapse response signal. 
4.3. Reservoir Identification, Correlation and Well-to-Seismic Ties 


The results for lithology and reservoir identification are presented in Figure 4. A total of nine sand bodies (A, B, C, 
D, E, F, G, H and I) were identified and correlated across all seven wells in the field. Two reservoir sands were 
selected for the purpose of this study (Reservoirs A and I). The resistivity logs which reveals the presence of 


hydrocarbons were used to identify the hydrocarbon bearing sands. 
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Figure 2. Reservoir Production Rate Versus Days _ Figure 3. Effect of Reservoir Pressure Decline on 
of Oil Production Production 
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Figure 5. Synthetic seismogram generation and 
well-to-seismic tie conducted for Fuba Field using 
Well-1 Checkshot 


Figure 4. Well section showing reservoir 
identified and correlated across Fuba Field 


Figure 7. 4D Effect (Diffference) Betweeen Base 
and Monitor Interpreted Faults Displayed on the 
Variance Time Slice 


Figure 6. Base, Monitor and Difference Seismic 
Section for Inline 8590 Interpreted 
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Figure 8. Third Order Polynomial Velocity 
model Utilized for Converting Reservoir Surfaces Figure 9. Reservoir Surface for Depth Surface A 
from Time to Depth 
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Figure 11. Depth Residual Maps Generated from 
Figure 10. Reservoir Surface for Depth Surface I Surfaces Converted Using the 3 Order Polynomial 
Velocity Function for Reservoir A 
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Figure 12. Depth Residual Maps Generated from Surfaces Converted Using the 3™ Order Polynomial 
Velocity Function for Reservoir I 


On Figure 4, the sands are coloured yellow while shales are grey in colour. The results for well-to-seismic tie 
conducted on Fuba field using density log, sonic log and checkshot of Well-1 is presented in Figure 5. A statistical 


wavelet (ISIS time) was used to give a near perfect match between the seismic and synthetic seismogram. 
4.4. Fault and Horizon Interpretation 


The results for the interpreted faults in Fuba field are presented in Figure 6 shows both synthetic and antithetic 
faults interpreted along seismic inline. Faults are more visible along the inline direction because this direction 
reveals the true dip position of geologic structures. The variance time slice was used to validate the interpreted 
faults as shown in figure 7. The areas in black circles indicate the 4D response. All interpreted faults are normal 
synthetic and antithetic faults. A total of thirty-six faults were interpreted across the entire base seismic data while 


forty faults were interpreted across the entire monitor seismic data. Of the interpreted faults, only F1 (synthetic 
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fault) and F16 (antithetic fault) faults are regional, running from the top to bottom across the field. Hence, these 
faults play significant roles in trap formation at the upper, middle and lower sections of the field. In Fuba Field, the 
number, lengths, dips and orientations of faults are not equal in both base and monitor seismic sections suggestive 
of fault reactivation from hydrocarbon production. The results for the interpreted seismic horizons (Horizons A 
and I) are also presented in Figure 6. On these horizons, the fault polygons were generated and eliminated. The 


horizons were used as inputs for the generation of reservoir time surfaces. 


4.5. Depth Surfaces 


The result of depth conversion residual analysis is presented in Table 1. The depth residual is the difference 
between the depth values of the well top from each well and the depth value from the depth converted reservoir 
surfaces. The depth residual analysis revealed that surfaces converted using the linear velocity function had the 
largest residuals ranging from -31.60 to +61.67 and from -50.58 to +40.84ft in reservoir A and reservoir I 
respectively. This is closely followed by the residual values obtained with the 2™ order polynomial function. The 
third order polynomial function shows the least residuals, ranging from -6.69 to +6.61 ft in reservoir A, and -9.48 
to +8.42ft respectively. The negative depth residual indicates that the depth conversion process displaces the 
reservoir to a greater depth than where it occurs in the subsurface, while a positive depth residual signifies that the 
depth converted result has placed the reservoir at a shallower depth (Ogbamikhumi and Aderibigbe, 2019). The 
resultant depth residual values generated using the various velocity models (linear, 2" and 3™ order polynomials) 
were compared in order to select the most suitable velocity model for depth conversion of the reservoir surfaces. 
Figure 8 shows the 3™ order polynomial velocity model which was selected and used as most suitable velocity 
model for converting A and I reservoirs from time to depth because it has the least residuals. The depth converted 
reservoir A and I surfaces are presented in Figures 9 and 10 for the third order polynomial velocity function. The 
depth structure maps reveal that the reservoirs are anticlinal and fault supported. Reservoir A is found at a 
shallower depth from 6500 to 7500 ft while reservoir I is found at a deeper depth ranging from 11500 to 13000ft 
respectively. The depth residuals recorded from the various well locations were used to generate depth residual 
maps which are presented in Figures 11 and 12 respectively. The depth residual maps revealed that higher 
residuals on reservoir A and I surfaces are associated with the eastern and western regions which are areas not 


penetrated by any well. 


Table 1. Depth Residual between Well Tops and Resultant Depth Surfaces 


Depth Depth Depth 
Difference Difference Difference 
‘ Surface Surface Surface 
Reservoir/ | Well Top (ft) (ft) (ft) 
Well (ft) (ft) (ft) (ft) 
Linear Velocity Function | 2nd Order Polynomial | 3rd Order Polynomial 
Reservoir A 
-7054.07 -7079.08 25.01 -7032.58 -21.49 -7053.13 -0.95 
Well-1 
Well-2 Missing Missing Missing Missing Missing Missing Missing 
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Well-3 -6877.06 -6849.08 -27.98 -6886.40 9.34 -6880.86 3.80 
Well-4 -6977.93 -7039.60 61.67 -7004.73 26.80 -6971.24 -6.69 
Well-5 -6905.39 -6873.79 -31.60 -6859.58 -45.81 -6900.65 -4.74 
Well-6 -7065.18 -7105.10 39.92 -7028.91 -36.27 -7070.87 5.69 
Well-7 -6846.24 -6877.44 31.20 -6854.65 8.41 -6852.85 6.61 
Reservoir I 

-11690.91 | -11720.12 29.21 -11674.22 -16.69 -11690.91 0.00 
Well-1 
Well-2 -11823.41 | -11780.54 -42.87 -11807.54 -15.87 -11823.41 0.00 
Well-3 -11650.06 | -11684.44 34.38 -11666.36 16.30 -11656.67 6.61 
Well-4 -11887.08 | -11845.26 -41.82 -11912.42 25.34 -11877.60 -9.48 
Well-5 -11599.86 | -11549.01 -50.85 -11581.29 -18.57 -11595.11 -4.75 
Well-6 -11569.00 | -11534.94 -34.06 -11586.60 17.60 -11564.27 -4.73 
Well-7 -11551.91 | -11592.75 40.84 -11534.64 -17.27 -11560.33 8.42 


== 5.0. Conclusion 


Production decline of 1192.21bbl/day resulted from pressure decline of 95.50bar/year or 0.062psi/year. Most of 
the data points were not repeatable as evidenced in the computed Normalized root mean square (NRMS) of 0.38 
meaning that on 62.0% of Base and Monitor data points were coincident. A total of nine sand bodies (A, B, C, D, 


E, F, G, H and I) were identified and correlated across all seven wells in the field. 


Two horizons (A and I) were selected for the study. Structural interpretation of seismic data revealed that the field 
is highly faulted with synthetic and antithetic faults which are in line with faults trends identified in the Niger 
Delta. Fault and horizon interpretation revealed that closures found are collapsed crestal structures bounded by two 
major faults. The depth structure maps reveal anticlinal faults. Reservoirs are found at a shallower depth from 6500 
to 7500 ft and at a deeper depth ranging from 11500 to 13000ft. The depth residual maps reveal higher residuals 


associated with the eastern and western regions which are areas not penetrated by any well. 


The synthetic and antithetic faults act as good traps for the hydrocarbon accumulation in the study area. The 
lengths, dips and orientations of the faults and horizons, in the base and monitor stacks, are not equal indicative of 
faults reactivation that could have resulted from hydrocarbon production. In reservoirs, hydrocarbons were 
encountered by all seven wells drilled in the field. Deterministic hydrocarbon volume estimation of the area of 
study indicates large quantities of hydrocarbon at the intervals where these structures were mapped. Within the 
limits of the available data, it is recommended that further studies should include integration of stratigraphic data 


of all the wells. This will provide more reliable data for interpretation of the depositional environments. 
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